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Abstract - In recent years, studies on cyclic molecular nanomagnets have captivated the attention 
of researchers. These magnets are finite in size and contain very large spins. They are interest¬ 
ing because they possess macroscopic quantum tunneling of Neel vectors. For antiferromagnetic 
molecular nanomagnets with finite number of even-numbered coupled spins, tunneling involves 
two classical localized Neel ground states separated by a magnetic energy barrier. The question 
is: can such phenomena be observed in nanomagnets with odd number of magnetic ions? The 
answer is not directly obvious because cyclic chains with odd-numbered coupled spins are frus¬ 
trated as one cannot obtain a perfect Neel order. These frustrated spins can indeed be observed 
experimentally, so they are of interest. In this Letter, we theoretically investigate macroscopic 
quantum tunneling in these odd spin systems with arbitrary spins s, in the presence of a magnetic 
field applied along the plane of the magnet. In contrast to systems with an even-numbered coupled 
spins, the ground state of the cyclic odd-spin system contains a topological soliton due to spin 
frustration. Thus the classical ground state is 2A-fold degenerate as the soliton can be placed 
anywhere along the ring with total Sz = ±s. Small quantum fluctuations delocalize the soliton 
with a formation of an energy band. We obtain this energy band using degenerate perturbation 
theory at order 2s. We show that the soliton ground state is chiral for half-odd integer spins and 
non-chiral for integer spins. From the structure of the energy band we infer that as the value 
of the spin increases the inelastic polarized neutron-scattering intensity may increase or decrease 
depending on the strengths of the parameters of the Hamiltonian. 


Introduction -. Quantum antiferromagnetic spin 
hain is one of the subjects in quantum magnetism en- 
owed with rich and interesting physics Effil- In re- 
ent years, a new class of quantum antiferromagnetic spin 
hain has emerged in molecular wheels, which comprises 
;iolecules with an even number of large magnetic spins ar- 
* anged periodically on a ring [5]|^ . These molecular wheels 
have finite length a few magnetic ions (sites) coupled an- 
tiferromagnetically. This makes their behaviour different 
from the conventional quantum spin chain [^|^, which 
usually contains small magnetic spins and large number of 
spin sites (thermodynamic limit). The two remarkable dif- 
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ferences between molecular wheels and conventional quan¬ 
tum spin chain are as follows. (1). Molecular wheels pos¬ 
sess macroscopic quantum tunneling of the Neel vectors 
separated by a magnetic energy barrier. (2). Molecular 
wheels possess neither magnetic phase transition nor anti¬ 
ferromagnetic long-range order due to the shortness of the 
ring (chain) and the presence of very large spin length. 
However, as the length of the spin decreases and the sys¬ 
tem goes to the thermodynamic limit, these effects might 
be observed. Thus, these molecular wheels provide ade¬ 
quate knowledge of quantum and classical regimes. 

In the last few years several molecular nanomagnets 
have been synthesized and studied experimentally by in¬ 
elastic neutron scattering method j^[^. There have been 
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numerous applications of molecular magnets, however, 
molecular nanomagnets with odd-numbered coupled spins 
pose serious challenges and spin frustration effects might 
be observed j^. There are few odd rings; recently some 
odd-numbered magnetic rings have been synthesized |10| 
and studied by inelastic neutron scattering and exact 
numerical diagonalization |11| . This includes Cr 7 with 
s = 3/2 and N = 7, which has open boundary condition 
and thus there is no spin frustration. Numerous studies on 
molecular nanomagnets have focused on the pure isotropic 
Heisenberg model |8PPI . This is usually a first approxi¬ 


mation because most magnets have dominant exchange in¬ 
teraction. The inclusion of anisotropy and magnetic fields 
can induce exotic phenomena. 

In this Letter, we consider the general Hamiltonian of 
a cyclic molecular magnet, modeled as N arbitrary large 
spins, regularly spaced on a circular ring [^; 


N 


N 


N 


H = Sz+1 -DY, Si + 9^iBB -Y^i- (1) 


1=1 


1=1 


1=1 


We assume periodic boundary condition S^r+i = Si and 
the case of odd N. The first term represents the mag¬ 
netic exchange coupling between spins with J > 0, the 
second term is the easy-axis magnetic anisotropy with 
D > 0, and the last term is the Zeeman coupling. This 
model has a wide range of applications in many cases of 
physical interests. For bipartite molecular wheels with 
even-numbered coupled spins, eq.Q is effectively a two- 
spin Hamiltonian |13| , which describes different molecular 
nanomagnets I® | |14| depending on the values of J and D. 
In the classical terminology, one describes such systems 
using spin coherent state (SCS) path integral formalism 
, where the spins are assumed to be confined on a two- 
sphere S 2 - In this formalism, the spins can be thought 
of as classical vectors aligned along the z-axis, but cant 
towards the direction of the magnetic field. The classi¬ 
cal ground states consist of two Neel vectors rii and n 2 , 
which are degenerate (for B = 0) separated by a magnetic 
energy barrier of height 2DNs‘^; see Fig.(0. When the 
barrier height is larger than the ground state energy, the 
Neel vectors are well localized at the two degenerate min¬ 
ima of the magnetic potential; quantum tunneling lifts the 
degeneracy of the ground states and produces an energy 
splitting, which oscillates in the presence of a hard-axis 
magnetic field (^[I^ {i.e. B = Bex). This formalism has 
been successfully employed to describe quantum phenom¬ 
ena in physical molecules |5|19| such as CsFeg, with N = 8, 


s = 5/2; NaFee, with N = 6, s = 5/2, etc. An exact nu¬ 
merical diagonalization also corroborates the results 
of SCS path integral formalism. However, for molecular 
wheels with odd number of magnetic ions, SCS path in¬ 
tegral is infeasible and the spin chain does not necessarily 
possess an effective two-spin Hamiltonian as the system 
is non-bipartite. In this case the classical ground state 
is frustrated and contains a topological soliton or domain 
wall, which is 2N-fold degenerate as shown in Fig. (©• The 


Wfold degenerate one-soliton states with total Sz = s can 
be written quantum-mechanically as 




ni=f =1 

Fig. 1: Schematic sketch of the two classical degenerate Neel 
ground states for antiferromagnetic molecular nanomagnets 
with even-numbered coupled spins in zero magnetic field, sep¬ 
arated by a magnetic energy barrier of height 2DNs^. 




Fig. 2: Schematic sketch of the 2A-fold classical degenerate 
ground states for antiferromagnetic molecular nanomagnets 
with odd-numbered coupled spins in zero magnetic field. 


\mu) = - ,t4, /M/ ,t4); (2) 

771,m+1 

where m = I,-- - ,A^; t= s, 4,= —s, and the subscript 
u denotes the up-up soliton. The corresponding Wfold 
degenerate one-soliton states \md) with total Sz = —s are 
obtained by rr-rotation of all the spins in eq.Q. 

In principle, molecular wheels with even-numbered cou¬ 
pled spins should contain topological solitons as excited 
states, but they must come in pair by hipping an odd num¬ 
ber of neighbouring spins from the fully anti-aligned Neel 
states. This has been demonstrated in molecular clusters 
with a pure isotropic Heisenberg model via an exact nu¬ 
merical diagonalization . We will consider a hard-axis 
magnetic field along the plane of the magnet and the limit 
in which the exchange interaction is small compared to the 
anisotropy term, i.e, {D 3> J). In this perturbative limit, 
we expect small quantum fluctuations stemming from the 
interaction and the magnetic field terms to lift these de¬ 
generate states; this results in delocalization of the soliton 
with a formation of an energy band. The non-perturbative 
limit {D ^ J) is tricky and requires thorough investi¬ 
gation for finite odd-numbered rings as the semiclassical 
approach is infeasible. For infinite chains, this limit corre¬ 
sponds to the Haldane limit . Recently, there have been 
some studies on short odd-numbered magnetic rings with 
pure Heisenberg exchange interaction and open boundary 
condition [^. To the best of our knowledge, no solution 
has been reported for short odd-numbered magnetic rings 
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with periodic boundary condition in this non-perturbative 
limit. 

Perturbation theory-. Quantum dynamics in short 
antiferromagnetic molecular wheels (chains) are usually 
studied by semiclassical formalisms and exact numerical 
diagonalization techniques. However, the basic principles 
of perturbation theory apply to both classical and quan¬ 
tum systems; it does not rely on the length of the spins 
(large or small) or the size of the systems (finite or in¬ 
finite). Thus, in many different occasions perturbation 
theory often reproduces the semiclassical results. We will 
follow the path of perturbation theory; it is customary to 
write eq. 0 as a sum of three terms: 


fluctuation next to the position of the soliton. The small 
fluctuations (splitting) from Vi at order 2s give: 

('S';;+i-S'++2)^'’ \mn) OC |m„ -f 2); 

('S';J;_iS'™)^" \mu) OC |to„ - 2); (10) 

where the states \mu ± 2) are defined periodically, thus 
they are elements of |to„). However, fluctuations from V 2 
at order 2s map to a state with total Sz = —s: 

\mu) OC \md - 1 ) ; |to^) oc \md + 1 ) ; 

( 11 ) 

where 


H — Hq + Vi -I- V2- (3) 

These terms are defined as 

N 

Ho = -DY,SI-, (4) 

^ / 1 \ 

+ 2 (5+^r+i + j; (5) 

h ^ 

^2 = Y E + E) ; h, = (6) 

where ^ raising and lowering spin 

operators. The TV-fold classical degenerate ground states 
are also degenerate ground states of Hq: 




m—l,m 

( 12 ) 

It is evident that fluctuations away from the soliton, such 

,S-,V^^nd(S+S- ^ 2 . 


i^m+4^m+5) 


and (-S')))_|_ 2 S'^_|_ 3 )^®, produce states with 


higher energy, which give small corrections in the pertur¬ 
bative series, thus can be ignored. 

Degenerate perturbation theory has been developed for 
single molecule magnets 22 . We employ this formal¬ 


ism to our problem; the matrix elements for the present 
problem are given by 


, (13) 

= {md\%'Bl^~^\mu), (14) 


Ho \mu) = So \mu ); (7) 

where £0 = —DNs^. At first order in degenerate pertur¬ 
bation theory, we must diagonalize V in the degenerate 
subspace |to„): 

•^m' ,mu “ {Bn^\V\TTlu) = Vn'^,m^ I ~ I; 2, • • • , N. 

( 8)1 


The splitting terms (raising and lowering operators) in 
V = Vi + V 2 map |m„) to states out of the degenerate 
subspace; consequently, the overlap between the resulting 
states and the degenerate states |m(j) vanishes. Thus, only 
the Sz term in Vi gives a nonzero energy contribution and 
the eigenvalues are determined by the secular equation 
(2§: 

<iet{Vrn',m “ = 0, (9) 

where Vm',m = Js^ {1 — {N — 1)) 6m',m', hence, the deter¬ 
minant in eq.([^ is indeed diagonal, yielding = S^ + 
£1 = -DNs'^-Js^{N-2). Evidently, A = = 0. 

A nonzero energy splitting is obtained when the transverse 
operators map the degenerate subspace into itself. This 
occurs at order 2s; hence, A = £;C«<2s) _ g{n< 2 s) _ 
However, at order 2s, we must consider small quantum 


where 

is the (2s — l)-th order transition matrix and Vc = ^ — 
X)m=i 1”^) (’^l complementary projection operator. 

The total matrix elements can be written as 


'vty? ‘2,s _ 0/^ 7 

— V/Crn, 




(16) 


Energy band-. The diagonalization of the Nx N matrix 
in the preceding section gives the energy band (spliting) 
of the system. The components of this matrix can in¬ 
deed be computed exactly for the one-soliton states. The 
only terms that survive in eqs.( 13 \ and (14 \ are indeed the 
fluctuations next to the soliton, as defined in eqs.(10l and 
dTTl). Thus, 


— (^lil (^^m-|-l,m+2 ^Trn — l,m ) |to„) ; (17) 

= {rh'dl {Tm + Tm+i) |m„) , (18) 


where the T operators are defined in the Appendix. The 
matrix elements are found to be: 




, 771 ^— 2 ] 

^h,j. [6rh'^,rhci + l 6m'^,rhfi — l]i 


(19) 
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where the are given by 




2s 2s 


n a(2s - fJ + 1) TT 

2Da(2s — a) 

(7 — 1 ( 7—1 ' ' 


= ASZJs^ 


J 

Id 

2s 2s 


2s 


( 20 ) 

( 21 ) 


=± ( ^ ) n \/f^(2s - o-+1) n 


a—l 


= ± 


2Ds 
(2s- 1)! {2D 


±DJ- 

TT 


Da{2s — a) ’ 

( 22 ) 

- / 7 \ 2s 

s I ehx \ 


^Ds ) 


■ (23) 


The last result in eq.(23l follows from the Stirling ap¬ 
proximation n\ ~ \/2Tm (") for s 3> 1. The first two 


products in eqs.(20) and (22) emanate from the relation 


of the raising and lowering operators. The second prod¬ 
ucts stem from the energy denominators. There is a plus 
or minus sign in these results which can be understood as 
follows. There are 2s — 1 negative energy denominators in 
eqs. and ([T^. Thus, if s is integer, 2s — 1 is odd, and 
a minus sign appears; whereas for half-odd integer spins, 
2s — 1 is even, and a plus sign appears. Due to the cyclic 
symmetry of the ring, in matrix representation we obtain 
a circulant matrix |24|: 


■ 0 



0 



"^/ix 

0 

^/tx 

••• 

0 




0 

'^/tx 


0 



^/tx 




0 

0 









0 


X/lx 


0 

... 


0 


_ 


(24) 

The eigenvalues are given by the finite cyclic group 
Fourier transform (241: 


7V-1 

em=Y. 

j=0 


C ■ 


^3 ^ 

™ y/N 


/47rm\ /27rm\ 

= 2^j cos () + 2%, cos () ■ 


It should be noted that the result of degenerate perturba¬ 
tion theory is unchanged if one uses the normalized linear 
superposition of the 2iV-fold degenerate one-soliton states 
|ws) = ^ (|wu) + \md)) with = 0. 

Chiral States and Polarized Neutron Scattering. 
In terms of the momentum q = 2Trm/aN, the energy 
bands can be written as 

Cq = —2cos(27rs)[|‘^j| cos (2qa) + \ cos (ga)]. (28) 

For spin-4 systems such as CsCoBra and CsCoCls with 
odd number of sites, the energy band reduces to a well- 
known Villain result [T6p8 . This is expected as large spin 
systems should in principle reproduce the results for small 
spin systems at least in perturbation theory. The major 
difference, however, is the size of the systems, N. Con¬ 
ventional quantum spin systems usually have large num¬ 
ber of sites (thermodynamic limit) and strong quantum 
fluctuations as opposed to molecular nanomagnets. As a 
consequence, spin systems with finite N and large spins 
usually behave as classical systems but possess quantum 
phenomena [T^. The one-soliton energy band in eq. (281 is 


thus the general formula for any system size and arbitrary 
spins. The argument qa and 2qa signify the total hopping 
sites of the soliton, as a result of small quantum fluctu¬ 
ations. For a single large spin, J = 0, eq.([^ reduces to 
H = —DS^ + hxSx, which describes Mni 2 Ac; then eq.( [2^ 
becomes the exact ground state tunneling splitting |21| 
with q = 0, which corresponds to the transition |t) <->■ |i). 
In order to characterize states, it is expedient to introduce 
the vector chiral operator defined as K = S; x Sj+i. 
The transverse component can be diagonalized in the one- 
soliton states |ms), where \ms) is the linear superposition 
of the 2N-fold degenerate one-soliton states. We find that 
K.^ = {rh'^\{KxI^\ms) is given by 


IC^ = sin(ga), 


(29) 


(25) 


where m = 0,1, 2, • • • ,N — 1 modulo N 

is the shift quantum number [^, and Cm are the row or 
the column entries of . The corresponding normalized 
eigenvector are column vectors with components 


where = 2(2s)!(s)^®. As shown in Fig.(|^, for half-odd 
integer spins, the soliton ground state is chiral located at 
g = ±^, while for integer spins the state at g = 0 is the 
non-chiral soliton ground state. 

Due to chirality the quantity of interest is the polarized 
neutron scattering intensity, which is an indispensable pro¬ 
cess in experiments |17| . For spin 1/2 systems the emer¬ 
gence of soliton in quantum antiferromagnetic spin chain 
has been studied experimentally using polarized neutron 
scattering techniques 18 . In the present problem, experi- 


(26) 


The one-soliton energy band for our model then simpli¬ 
fies to 


mental measurement should be reminiscent of that of spin 
1/2 systems [^. The crucial difference, however, is that 
in the present problem the energy band is formed at order 
2s and the system size is finite and short. The quantity of 
interest is usually the magnetic scattering intensity, which 
has the form 18 : 


(27) 


^ 4 ± 
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Fig. 3: The one-soliton energy band for different values of 
the spins with JjAD — 0.3 and hx/2D — 0.05. The ground 
state for half-odd integer spins is chiral located at g = ±7r/2a; 
whereas that of integer spins is non-chiral located at g = 0. 
As the spin increases the energy band approaches a constant 
value. 


soliton. As the soliton can be placed anywhere along the 
cyclic chain, the resulting energy is A^-fold degenerate in 
the sector with a total Sz = s. Small quantum fluctuation 
stemming from the interaction term translates the soliton 
by two magnetic sites in the same state, whereas fluctua¬ 
tion from the magnetic held term translates the soliton by 
one magnetic site in a different state with total Sz = —s. 
Delocalization of the soliton by quantum fluctuation forms 
an energy band at order 2s in degenerate perturbation 
theory. We obtained the energy band and showed that 
the soliton ground state is chiral for half-odd integer spins 
and non-chiral for integer spins. We further demonstrated 
the structure of the inelastic polarized neutron scattering 
intensity, which should be of experimental interest. Our 
results are general and apply to systems with arbitrary 
size and spins. 


* * * 


where Q = q — q' is the momentum transfer, oj = E — E' 
is the energy transfer, and a,l3,'y = x, y, z. The plus or 
minus sign denotes the directions of the magnetic field. 
The magnetic dynamical structure factor S°‘^{Q,u}) and 
the dynamical chirality /C'’'(Q,a;) are given by 


S°‘P = 


= 




Q 




S°‘^. 


(31) 

(32) 


where S'q“ = •S'q — Qa is the molecular spin 

component transverse to the momentum transfer and Cap-y 
is the antisymmetric Levi-Civita tensor. In the absence 


of chirality, eq. (301 reduces to usual formula for unpolar¬ 


ized neutrons scattering |6[|25|. The thermal average in 


eq.(31l must be performed in the one-soliton states |g). 


which represents coherent superpositions of solitons at all 


the magnetic spin sites. As shown in eqs.(lO) and (111, 


a non-zero correlation function is obtained by taking the 


transverse components in eqs.(31) and (321 to order 2s. 


This is the case for every quantity that requires the com¬ 
putation of correlation functions such as the spin-lattice 
relaxation time. Without computing the neutron scat¬ 
tering intensity explicitly, we can get an intuitive under¬ 
standing of its structure. Depending on the strengths of 
the magnetic field and the coupling constants J and D, 
the energy band tq may increase or decrease as s increases 
(see Fig. @). Thus, it is evident that the scattering inten¬ 
sity will have a similar trend but approaches a constant 
value at very large s. 

Conclusion—. In conclusion, we have studied an 
antiferromagnetic molecular nanomagnet with an odd- 
numbered coupled spins, in the presence of a transverse 
magnetic field applied along the hard-axis anisotropy. We 
considered arbitrary spin with periodic boundary condi¬ 
tion. In this case we showed that the Neel state is frus¬ 
trated, resulting in a defect which introduces a topological 
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Appendix-. The T operators in eqs.([T7| and (18) are 
given by 

2s / \ 2s—1 




m+l,m+2 


= (w) 






Tm = 


£o — Ho 
P, 


C- 0+ 
^m+l^m+2 


£o — Hq 


^0 — Ho 


-s- 


q-i- c- 


(33) 

1 

? 

(34) 

(35) 
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